Stability of the 3-form field during inflation 
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We consider the minimally coupled 3-form field which has been considered as a candidate to realize 
inflation. We have studied the conditions to avoid ghosts and Laplacian instabilities and found that 
some classes of potentials, e.g. the Mexican-hat one, will in general be unstable. We propose other 
classes of potentials which are instead free from any instability, drive a long-enough slow-roll regime 
followed by an oscillatory epoch, and as a consequence, can provide successful inflation. 



I. INTRODUCTION 



o 

(n: . . . . . 

, The inflation paradigm was introduced in 1980 as a way to solve different issues, namely: the magnetic monopoles, 

■ the flatness, and the horizon problems [1]; however, it can also account for the observed temperature anisotropics in 
[ the Cosmic Microwave Background (CMB) [2] as well as the galaxy power spectrum [3|]. In other words a sufficiently 

long stage of accelerated expansion has been proposed as a way to solve all these problems at the same time. In 
' order to explain this period of accelerated expansion, some new physics is introduced, and a scalar field P, |3| (or 
more than one [5]) is commonly used. However, the real mechanism for inflation is yet unknown, so it is interesting 
' to explore different possibilities which in general may lead to different predictions for several inflationary observables 
Q-i (i.e. spectral index, tensor-to-scalar ratio 0,1^, non-gaussianity parameter ■ 

I Since fundamental scalar fields have not been discovered yet in nature, the idea of inflation might well be realized 
1^ . by other, higher-form fields. For example, vector inflation (or one-form inflation) has been intensively investigated 

■ Unfortunately, most of the vector-field models encounter instabilities [17H2(J |. More in general, the N-forni 
field inflation has been also investigated (2ll . [2^ and one of the results is that one-form and two-form fields are not 

.1 stable, whereas the three-form fleld can be stable [l^. Note also, that the four-form field models correspond to the 
>; /(i?) theories |2l,[2|. 

\^ , Recently, a form of inflation based on the evolution of a 3-form has been studied (25l - l28| . The origin of such a 
OS- non-standard form for the inflaton may come from a high-energy-scale theory, such as string-theory. Indeed it is 
interesting to study such a model, as it may provide an alternative way to obtain inflation. Since the essence of 
the 3-form is by construction different from a single scalar field, we expect this difference to play some role both at 
OA I background and perturbation levels. 

■ In fact, in this paper we will study the stability of a minimally coupled 3-form during inflation with a general 
] expression for the potential. We will then find the conditions which avoid ghosts and Laplacian instabilities (i.e. 

. . . we require a positive kinetic term and a non-negative speed of propagation for the independent linear perturbation 
modes). Once these conditions are obtained, we reconsider some models which have been recently introduced [21, 2^, 
and show that, if the potential is not carefully chosen, both ghosts and Laplacian instabilities will occur. 

Hence, we provide some classes of potentials, which, by construction, are instead free from these instabilities, and, 
\ in this context, we study their background evolution, in order to confirm that a slow-roll period of inflation is then 
followed by a regime where the 3-form oscillates, ending inflation. We will discuss the details of reheating, and the 
bounds on the inflationary parameters (spectral index, tensor-to-scalar ratio, and non-gaussianities) in a future work. 

The paper is organized as follows. In section [Hi we introduce the Lagrangian of the model and write down the 
equations of motion. Linear perturbation theory for this model is studied in section IIIIl where we give the no-ghost 
conditions and the squared speed of propagation for the scalar, vector and tensor modes. We present some classes of 
potentials which make the model free from ghosts and Laplacian instabilities in section IIVI where we also show that 
a slow-roll period of inflation is followed by an oscillatory regime which ends inflation. We write our conclusions in 
section |Vl 

II. THE MODEL AND THE BACKGROUND EQUATIONS OF MOTION 

Let us start with the following action 



S — d Xy^—g 



2 48 



(1) 
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where Aa/s^ is a 3-form, and F = dA is its Maxwell tensor [2J| , whose components can be written as 

A. The background 

Let us now consider a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) manifold whose metric element can be 
written as 

ds^ = ~dt^ + a{tfdx^ , (3) 
and, on this background, considering Eq. (jAll) in Appendix, the background 3-form A^p-y can be written as 

Aoij = , Aijfc = a^eijfe X , (4) 
where tijk is the three-dimensional Levi-Civita symbol (with £123 = 1). Let us define the following quantities 

V = V{y). F,^^^, Vyy^^^^, where y ^ A^0^A"f'\ (5) 

On a FLRW background, we have y — 6X^. As a consequence of these definitions, on FLRW, we have 



V{-X) = V{X), (6) 

V = UXXVy (8) 

I d_ (V^\ ^ XV XX - Vx 

12X dX V 12X / 144X3 



so that we will restrict the form of the potential V to even functions oi X . In this case the Friedmann equation can 
be written as 

El = UI^iH^ -px^O, (10) 

where 

px ^^X^ + V + ^H^X^ + 3HXX = ^Y^ + V , (11) 

is the effective energy density of the 3-form, and we have defined Y = X + 3HX. The second Einstein equation reads 
as follows 

E2 = A4i{2H + 3H^)+Px^0, (12) 
where px is the 3-form effective pressure defined as 

Px = -(^\x' + V + 3HXX + ^ H^X^ - l2VyX^^ ^2Vyy-px. (13) 

The equation of motion for the field gives 

Ex=X + 3HX + 3XH + 12VyX = Y + 12X Vy = 0. (14) 
The equations of motion are not all independent, due to Bianchi identities: indeed we have 

Ei+3H{Ei- E2) + YEx (15) 
One consequence of the equations of motion is 

M^iH = -Vyy, (16) 
so that the universe will be super-accelerating when Vy < 0. 
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III. LINEAR PERTURBATION THEORY 



A. Scalar modes 



Let us consider now the metric for the scalar perturbations in the following form [2£ 

ds^ = -{1 + 2a)df + 2^^■^|J dt dx' + a^il + 2<^>)dx^ , 



(17) 



where we picked a spatial gauge so that the three-dimensional metric is diagonal. As for the 3-form, by using once 
more Eq. (|Al[l given in Appendix, we can use a time gauge to fix the scalar perturbations as [27| 

(18) 



Aoij = aeijkdkl3{t, x) , Aijk = a^e^fe X{t) . 
By expanding the action at second order in the fields we obtain 

5(2) , 



dtd'xa^l^^:^{dijf - 2A/|i(ffa - $)^ 

2 a4 
1 



+ {Ya + 12KyXV' + 3r$) 



QM^yH^ - 2M, 



PI 2~~ 



(y2 



PI 2~ 



3(r^ + i2y.yX2)$ 



12KyX^ 



(19) 



At a first look, this action has important differences with the general action (for the perturbations) of scalar tensor 
theories (s^l- First of all the presence of the terms (9-0)^ and which, for a second-order general scalar-tensor 
theory, vanish after using the equations of motion. Both these terms now vanish only when the 3-form is absent. 
Furthermore the field /3 is not dynamical and it can be integrated out in Fourier space (together with a and t/;). 

In order to remove these auxiliary fields it is convenient to work in Fourier space: in this case, we can integrate out 
the fields a, V'j and /3, by using their own equations of motion. In Fourier space, with $(t, x) = (27r)~^/^ J d^k^i^e^^'^ ^ 
with reality condition <^_f. = the equations of motion for the constraints give 



UVyX'^ip + 2M|i(iJa - $) - UVyXP = , 



(20) 



and 



2M2ji?fc2^ 



fc2<I> 

QMl^H^ + 2M^i— + 3(r^ + 12VyX'^)^ 



Yk^l3 



0, 



(21) 



■ /3 -I- 12Vy/3 -Ya- UVyX^P - 3^$ = , 



(22) 



where we omitted the tilde of the Fourier modes for simplicity. This last equation can be solved for /3 as 

a'^{Ya + 12VyX4> + SY^) 



so that we also have 



^'-'Pl 

6VyX^ 



f3- 



2a^M^i 



fc2 + UVyu'^ 



3a^Y 



6VyX^ 



a'^{2Ml^H - XY) 



and finally 



M^^HP + 6 M^^V^ya^{3 HX^ + 2 Ml^H - XY) 
~ Ml^H [Ml^BH + 6 V^yo? (3 HX"^ + 2 Af2^iJ - 2 XY)] 
6 VyMl^k'^X^ + 18 VyO^X (6 X^^y -I- Ml^HY) 
^ A/2ji/ [Ml^k^H + 6 yj,a2 (3 i/A:2 + 2 Af^ji? - 2 XY)] 

By substituting these expressions into the action written in Fourier space, we find 

• • A:2 

$fc$_fc - c^(t, fc^) $fc$-fc 



5= / dtd-'kQ{t,k^) 



(23) 



(24) 



(25) 



(26) 
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The no-ghost condition is then found to be 



where 



Ghost conditions 



>0, (27) 



It should be noticed that condition ([?7)) should hold at all times during inflation, whether or not the trajectory is in 
a slow-roll regime. On using the equations of motion (by replacing AIp^H^ with the Friedmann equation and then 
3HX = r - X) we find that 

3 H^X^ + 2 Af|iij2 - 2 XYH =^X^ + -V, (29) 



3 3 



so that 



^ 6aH4,V,Y^ 

M^.^k'^m + 2 Vya^ (I2 + 2V)' 

This quantity must be positive for all /c's. For high fc, we find the condition Vy > 0. This condition must be 
satisfied along the trajectory of motion. In some cases, it may be possible that for some (positive) values of y, Vy is 
negative, but such values of y are never reached: in this case the model can still be viable. We note here that the 
condition Vy > 0, on using Eq. (fT6|). forbids the dynamics to be super- accelerating. 

For low fc's we find another requirement, that is X^ + 2V > 0. Once more, this condition must be satisfied along 
the trajectory of motion. The bottom line is that the two conditions > 0, V^j^ > are sufficient conditions for not 
having ghosts. If these conditions are not satisfied for all (positive) y's, one should check that, at least for the values 
of y along the trajectory of motion for the model, the above mentioned conditions still hold. 



2. Speed of propagation 

The speed of propagation is found as the large-A: limit of c^(t, fc^) of Eq. ([26|) . One can show that the speed of 
propagation, on using the background equations of motion, is given as 

4= limc2,(t,fc2) = i + ^ = ^. (31) 

fc-i-OO V^y V^X 

The speed of propagation found here corresponds to the one found by Koivisto and Nunes In general, only the 
simple quadratic potential V on y^ implies a propagation with speed of light for all dynamics. Since y > 0, then a 
sufficient condition to avoid also Laplacian instabilities (besides the ghosts, Vy > 0) is Vyy > 0. 



B. Vector modes 



Let us define the metric perturbation for the vector modes as 

6goi ^ aGi, and Sg^j = a^iCi^. 



(32) 



where Gi^i = = C^.i. We will also choose a gauge for which the 3- form has no vector perturbations (uniform field 
vector-gauge). This choice completely fixes the gauge degrees of freedom. In this case, one can show that the action 
for the vector modes becomes 



I 



dtcfs 



Ga^VyX^C^Ci + Ua'^VyX^C^Zi + ^M^-^a {djZ,){djZ,) + 6a^VyX^Z,Z, 



(33) 
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where we introduced the field Zi = Gi — aCi. By introducing Fourier modes, it is possible to integrate out the field 

Zjj cLS 



so that the action for the vector modes becomes 



5* = J dtd^kQvit,k^) d^it,k)d^{t,~k) , (35) 
so that it is clear that the vector modes do not propagate. 

1. Ghost condition 
The ghost condition for the vector modes corresponds to Qv > 0, that is 

^^-MpM"24^n^>°' ^^^^ 

implying 

Vy>0, (37) 

which coincides to one of the conditions already found for the scalar modes. 

C. Tensor modes 

The tensor modes are not affected by the presence of the 3-form, as this latter one is minimally coupled to gravity 
and it does not possess tensor degrees of freedom. To show this more in detail, we choose the tensor perturbations as 
6gij — hjj — ft-+e^ + hye^^^ where both the symmetric tensors e^j are transverse and traceless. We also impose the 
normalization condition, eij{k)eij{—k)* = 1, for each polarization, whereas efj{k)e^j{—k)* = 0. Therefore the second 
order action can be written as 



St= [ dtd^xa^ 



M2 



hi--idh,r 



(38) 



A=+,x ' 

SO that no stability condition comes from the tensor sector. 

IV. SUITABLE FORM OF POTENTIALS FOR 3-FORM INFLATION 

According to the previous section, one of the no-ghost condition can be written as V^y > 0, where y — 6X^ > 0. 
The existence of ghosts in the model depends on the shape of three form potential, but not on the sign of X (as 
y oc X^). In fact, in order to search for the form of potentials, which makes the 3-form field ghost-free and without 
Laplacian instabilities (cx > 0), we need to study more in detail the evolution of y (or, equivalently, X). From the 
Friedmann equation, we have 

so that the 3-form field can play the role of a slow- rolling inflaton if Vx X / Mp^ ^ H^. Substituting the above Eq. ((39)) 
into the evolution equation (fT4l) . we get 

X + 3HX + Kff.jf = , (40) 

where 

dV^s A, 3 X2 
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so that the effective potential is given by 



3 e 



2 



On using the dimensionless variables 



X , 3.T + x' 

X = , and w = 



where a prime denotes a derivative with respect to = In a, Eq. (j40p can be written in the autonomous form as 



a;' = 3 



-w-x 



where we have introduced the function 



w^' = ^ (1 - ixw-J'^ 



In these variables the slow roll parameter can be written as 



H 



X{l-w^)x. 



(42) 
(43) 

(44) 
(45) 

(46) 

(47) 



The accelerating expansion of the universe is acquired by demanding e ^ 1. From this parameter, one can see that 
the kinetic term does not necessarily need to be small compared to the potential term, as for the standard picture 
of the inflaton scalar field. Conversely, it requires that « 1 when x X{x) ^ 0{1). In order to have inflation, one 
needs one more requirement to guarantee that the accelerating expansion is long enough. We introduce a parameter 
to characterize this behavior as 



l-2e 



(1 



(48) 



where the inflationary period requires that \ri\ <C 1. Since > 0, \ri\ will be small if x' /x is small, that is x needs to 



be in a slow-roll regime. From Eq. (|44p. it implies that jryj <C 1 will be satisfled if x ~ y 
have imposed also the flrst slow-roll condition |e| <C 1 when x X{x) ~ 0(1). 

We note that, as a consequence of the deflnition of w, we have Y/{MpiH) — \/6w. 
equation (jlOp implies 



1 = px/{iM^iH^ 



w 



2/3ti; ~ ±-\/2/3, where we 
Therefore, the Friedmann 

(49) 



so that, if > 0, then < < 1. 

As we have already said, the field will slow- roll when it reaches the points P = {x,w) — (±-^2/3, ±1) in phase-space 
because these points are (de Sitter) fixed points (unless X{x) is not finite at these points). There might be other fixed 
points, M, which correspond to the points where A vanishes, that is AI = {x,w) — {x, •\/3/2a;), where X{x — x) ~ 0. 
It can be seen from Eq. (PT|) that the points P and M are the values of X which correspond to the extrema of the 
effective potential. 



A. Stability of the fixed points 

Let us start by studying the stability of the fixed point P. By choosing x = ±^2/3 + Sx^^\ and w = ±1 + Sw^-^^ 
we can linearize the equations of motion with respect to the small quantities 6x and 6w, and we find 

d{5x^^'^)/dN = ^/&5w^^^ - 3fe(i) , (50) 
d{5w^^^)/dN = 0, (51) 
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with solutions Sw'^^'^ ^ bi = constant, fa^^^ = ^2/36i +626 and 61,2 are (small) initial conditions. The fact that 
cP{Sw^^'')/dN'^ = 0, implies that, at linear order, we cannot deduce whether the fixed point is stable or not. In order 
to check the stability of this fixed point, one needs to study also the second-order solution. Therefore, by re-iterating 
the method, and choosing x = ±-\/2/3 + Sx^^^^ + Sx^^\ and w = ±1 + Sw'-^^ + Sw'-'^\ we find 



d{6x^^'^)/dN = V6J«;(2) _ 3^2;(2) , 



which can be solved as 

5x^^^ = 7273^3 + A(±V273) 6^(4/3 - 4iV) + ^/l\{±y/2j?,) 6162 iVe'^^ -t- 646" 



3N 



(52) 
(53) 

(54) 
(55) 



where 63^4 — 0{h\). These second order solutions show that, at second order, the fixed point is unstable (due to the 
term proportional to N). The fact that this instability appears at second order means that the instability will in 
general evolve slowly. This slow instability will make inflation end eventually. 

As for the fixed point M — {x, w) = {x, y/3/2x), where \{x) — 0, then by choosing x — x + Sx, and w = y^3/2x+Sw, 
we find the linearized equations 



where 



r = 



5x 
Sw' 

'V 



V6 



V 



3fe , 
(2 - ixp- fVSx, 



V 



The solution leads to 



where 



fe = die-^(3+7)/2^^^g-JV(3-7)/2^ 



7= a/9 -3r (3^2 -2)2 



(56) 
(57) 

(58) 

(59) 
(60) 



An instability will appear if 7 > 3, or F < 0. 



B. Mexican-hat potential 



Let us start with the Mexican-hat type potential of the form 

V^V,{x^-c^f 



This potential yields 



v.. 



V 



X 



v.. 



12X 



12M2jX 



3M2j 



(a;2 - c2) 



and 



3a;2 



-X 



(61) 



(62) 



Hence, the ghost will not exist if > |c| and for this range of x the 3-form is stable, i.e., > 0. For this form 
of the potential, the effective potential has 5 extremum points at a; = 0, a; = ±|c| and x — ±-^2/3. If |c| < ■\/2/3, 
the points x — ±|c| are the the minimum while x — ±-\/2/3 are the maximum points of the potential. The points 
x — ±|c| and x — ±-\/2/3 become the maximum and minimum points respectively when |c| > •\/2/3. In Figure ([T]) 
the potentials for the case |c| < ■y/2/3 and |c| > a/2/3 are plotted. In the case where |c| < a/2/3, the field x that 
starts at points A and C usually evolves to point B. Roughly speaking, the field can climb up the effective potential 
from the point C to A because x' is negative (when the field is on the right side of the point A), and |w| is never 
larger than unity (as V in this range is always positive). 
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Figure 1. The potential V = Vo{x^ — (?Y is represented by solid line while the effective potential is represented by dashed line. 
In the left panel c= 1/2 < -\/2/3 while c = \/3/2 > in the right panel. 





Cy 



Figure 2. The evolution of x for the potential V = V(,{x'^ — c?)'^ , and c — 1/2. We chose initial conditions a;(0) ~ 0.816188, and 
w(0) ~ 0.99962, so that e(0) ~ 9.289 x 10"'*. We notice that as x becomes less than one half, then a ghost appears, and, at the 
same time, c?x diverges and changes sign. 



By using the previous qualitative analysis, we can generalize this behavior to other potential forms. Since e oc 
\{x) a; (1 — lip'), the field can drive inflation if \w\ ~ I and A(x) x is not large. From equation pS)) . the inflation will 
be long enough if \x\ > y^2/3 initially. The autonomous equations imply that x' < if x > -\/2/3 and a;' > if 
X < —a/2/3. When we initially put the field in the region |a;| > 1/2/8, the field will evolve to the de Sitter point P , 
where the condition \ri\ <C f is satisfied. However, since P is not a stable fixed-point, the field will evolve by crossing 
P, and entering the region |x| < ■y/2/3. If there is a minimum in the region \x\ < the field will eventually 

oscillate around this minimum, so that a mechanism to end the inflation is allowed. We note that condition \r]\ <^ 1 
can be satisfied even if we initially put the field in the region x < •\/2/3. However, we may need to fine tune the 
initial value of w close to unity. 

For this potential form, when the field reaches the point B it will oscillate around this point. Nevertheless, the field 
becomes ghost when it moves to the left of the point B, so that the field cannot oscillate around the minimum of the 
effective potential without becoming ghost. 

In Figure [5] we plot the evolution of x for c = f /2 and some initial conditions compatible with slow-roll, to explicitly 
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V(x) 




Figure 3. The potential V{x) oc is shown here. In the figure, the bare potential and the effective potential are represented 
by a solid and dashed line respectively. 



show, also numerically, that the ghost (and the Laplacian instabilities) will, in general, appear. 

For the case |c| > y^2/3, the right panel of Figure ([I]) shows that the ghost will not exist and > if the field x 
evolves between A and B. Nevertheless, the point F is the stable fixed point in this case, so that the field that starts 
from point B will evolve towards point F. This implies that, also in this case, the field will become a ghost, eventually. 



C. Power- law potential 

We would like to present here a class of potentials which can give a stable inflation scenario. On choosing a power 
law potential of the form V (x y^, we immediately notice that Q cx y^~^, which in general vanishes (for p > 1) or 
diverges (for p < 1) as y ^> 0, unless p = 1. Since we will focus on values p > 1, most of the potential will allow the 
field to cross this value {y — 0), so that Q will vanish in the origin. This property represents a problem, in general, 
as this means that, at that point, the second order Lagrangian vanishes (as Cx remains finite for V cx j/^), and the 
theory becomes strongly coupled, i.e. higher order corrections become dominant. In order to avoid this possible issue 
we propose the following generalized power law potential. 



y(x2) = Fo[(x2)P + 6a;2 



(63) 



where p is a constant which can be, as for now, positive or negative, whereas & > 0. For this form of the potential, we 
have that the no-ghost condition 



X 



2V^\p{x') 



fe] >0, 



(64) 



is always positive for p > 0, and also finite for p > 1. For this reason, from now on. we will consider only the case 
p > 1. On the other hand, since 



(2p- l)p(x^)P-i +6 
p(a;2)P-i + h 



(65) 



will be always positive and finite for p > 1. The bottom line is that for p > 1 the model is free from instabilities 
for any real value of x, that is for any dynamics. It should be noticed that for the value p = 1, the potential reduces 
to a quadratic power law potential. On the de Sitter fixed point P, we have 



c\(x 



(2p-l)p(2/3) 



p(2/3)P-i + 6 



> 1 



(66) 



^ In general, we can choose a larger class of potentials given as V{y) = Vb (cy + Ci y^^), where c > 0, Ci > 0, and pi > 1. 
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Figure 4. The evolution of x and e for the potential V = Vb(l + b) , and p = 1. We chose initial conditions w{0) = 0.99, and 
2;(0) = 3-\/2/3 w(0). When the field a; start to oscillate around the minimum, the parameter e will start to oscillate around 
e = 3/2. Thus the inflation will end when e ~ 1 corresponding to ~ 72. 



and the inequality holds for p > 1, and h> Q. 
Because of Eq. , it can be shown that 



.2p-l 



bx 



and 



Kff = Vo 



2(P+1) 



-[2+p(2-3x^)] + 6a;^ ( 1-- 



so that the extremum points of this potential occur at a; = ±^2/3, and a; = 0. We also notice that 



\{x) X 



2 [p (a;2)P-i + h\ 



\{x- 



2^p-l 



6] 



(67) 



(68) 



(69) 



which, for p > 1, is positive and finite for all x. Furthermore lim^j^o Aa; = 2. 

For illustration, we consider here the simplest case p = 1 whereas, in the next section, we will describe in more 
detail the case p = 2 (quartic potential) . The potential and effective potential for the p = 2 case are plotted in Figure 
[31 Similar to the case of Mexican-hat type potential, the field that starts at |a;| > ■\/2/3 with \w\ ^ 1 is able to drive 
long-enough infiation. However, this time, the field rolls down the potential from the points A or B and then oscillates 
around the minimum of the potential without ghosts (y,y = Vb(l + 6)/6 > 0) or Laplacian instabilities (c^ = 1). We 
also use a direct numerical integration to confirm both the slow roll and the oscillatory regimes as shown in Figure |31 
From the evolution of e in the right panel, inflation ends at ^ 72 corresponding to e ~ 1. 



D. Quartic potential 

We study here a particular case of the power law potential introduced in the previous section, namely 

1/ = Vo(a;'' + 6a;2), 



(70) 



and we plot it (together with its effective potential) in Figure [S] 
In this case, we also have 



\ — up- 



(71) 
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Figure 5. The potential V = Vb(a;* + hx'^), with h — 1 (the continuous and black curve), together with the eilective potential 
Vcfi (represented by the dashed curve). 
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Figure 6. The evolution of x and w for the potential V = Vo{x'^ + bx^), and h — 1. We chose initial conditions to(0) ~ 0.9991, 
and x{Q) ^ 0.8158. 



so that, by introducing a dimensionless cosmic time t, we can write 



dx MpiH dx 



bx^ dx 



dN 



(72) 



In Figure [SI we show the evolution for both x, and w. The field slow-rolls until an oscillatory regime starts, making 
inflation end. 

During the slow- roll regime the propagation speed takes the value w '^i+sb'* ' whereas, as the solution starts 
oscillating, — ^ 1. This behavior is confirmed in Figure [7l 

Finally, we show in Figure [5] that after inflation ends, there is an oscillatory regime which mimics a dust dominated 
universe as we have cx a^^ cx e^^^. This behavior is similar to the standard single- field infiationary models, and 
this is not surprising, because as a; — 0, we find Vcs « bx^; so that the equation of motion for the field, Eq. ([T4)) . 
reduces to 



X + 3Hx ; 



-bx . 



for |a;| < 1 , 



(73) 
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Figure 7. The evolution of for tlie potential 1/ = Vo(a;'* + ba::^), and h —\. The model does not possess any instability. 
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Figure 8. The evolution of 3Mpi_H"^/Vb for the potential V = Vo(a;* + 6a;^), and 6=1 (continuous black curve). This figure 
shows that after inflation ends (around N ~ 70), the universe enters a matter dominated epoch, as the curve approaches a 
dashed line, which represents the line ln(3Afpii/^/Vo) — —SN + constant. This means that after inflation oc e""^'^ oc a~^. 



which exactly matches the equation of motion for standard inflation in the presence of a quadratic inflaton potential. 
In other words, the dynamics of the 3-form, for a; — 0, tends to be more and more identical to the dynamics of a 
single scalar field oscillating around the minimum of a quadratic potential. 

After inflation ends, during the oscillatory regime, in Figure [HI we see that a; — > 0, whereas w oscillates between — 1 
and 1. Furthermore, we also find that dx/dt — )■ together with x, whereas dw/dt keeps oscillating, remaining finite, 
as shown in Figure IHl 



E. Gaussian potential 



We now consider the exponential potential 

V = Voie'^y"' - 1) = yo(e''"' - 1) , (74) 
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Figure 9. The evolution of dx/dt as a function of x, and dw/dt as a function of w for the potential V = Vo{x'^ + bx^), and 
6=1 during the oscillatory regime. We chose initial conditions so that a.t t — 1, the values of x and w correspond respectively 
to x{N = 72) « 0.0644, and w{N = 72) « -0.46 of Figured We stop the integration at t = 100. This figure shows that x 
spiralizes, whereas w keeps on oscillating during the matter dominated regime. 

where z/ is a constant parameter which can be positive or negative. For this form of potential, we have 

^ = 2,.Vo e^"^"^ , (75) 

X 

so that the ghost will not exist if v is positive. The speed of propagation in this case is given by 

c| = 1 + 2iyx'^ . (76) 
This implies that if the ghost does not exist, is always positive. Substituting Eq. ([74l) into Eq. (j4T|) . one gets 

Kff.. = 2jyxVoe''^" (^1 - ^ x^^ , (77) 

or 

Kff - - ^ {3 + + e''^" [(3x2 _ 2)i/ - 3] } , (78) 

and we plot it, together with the bare potential, in Figure [TOl We also notice that as 2: — ^ 0, then VcS — Vqi^x'^, so 
that we expect an oscillatory regime to take place, ending inflation. 

It is easy to see that the effective potential has the extremum at x = ±-\/2/3 and x = 0. Similar to the analysis 
for the previous potentials, the field can drive inflation when we initially put it in the region satisfying the condition 
e^l,e.g. |x|> •y/2/3 and \w\ ^ 1. The condition jvyj ^ 1 will be satisfied when the field is frozen nearly x — ±y^2/3. 
Since x = ±-^2/3 are not stable fixed points, the field can continuously evolve through x = ±y/2/3 and then oscillates 
about X = eventually. This behavior is also shown by using numerical integration methods as seen in Figure [TT] 
Because of this behavior, the speed of propagation will be approximately equal to w 1 + 4i//3 in the slow-roll 
regime, whereas — > 1, as a; — > 0. In Figure \V2\ we also show the behavior of the Hubble parameter during the 
oscillatory regime, confirming that a matter-dominated era takes place during this epoch. 

Finally, we show the trajectory of dx/dt and a;, together with dw/dt and w in Figure [T3l 

V. GENERAL CONSIDERATIONS AND CONCLUSIONS 

From the investigation of the previous section, one can see that the viable 3-form models can be characterized by 
the shape of their potential. The study of power-law potentials suggests that the viable potential form which is free 



14 




1.0 -0.5 0.0 0.5 1.0 



Figure 10. The potential V{x) = Vo (e"^ — 1) is represented by a solid line (for u = 1), whereas the effective potential is 
represented by a dashed line. 
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Figure 11. The evolution for x (left panel) and e (right panel) for the potential Vo{e"^ — 1), and i/ — 1. We chose initial 
conditions w(0) ^ 0.999, and x(0) « 0.815. 



from ghosts and Laplacian instability should have the local minimum at x = 0. This is because when changes 
sign around the minimum point (as V is, by construction, an even function of x, V{x) = V{—x)), x also changes sign 
such that Vx/x = 12Vy is always positive (where y = 6X^, and x = X/Mpi). In this situation, cj^ > around x = 
because Vxx > 0. The speed of propagation is still positive as long as x remains significantly far from the nearest 
local maximum (if it exists) of the potential along the trajectory of motion. Hence, if the bare potential has no local 
maxima between x = ±-y/2/3, the field can evolve between x = ±-\/2/3 (which, in turn, are the local maxima of the 
effective potential) without giving rise to ghosts or Laplacian instabilities. 

In section HVl we have proposed a class of potentials which are free of instabilities, can drive inflation, and provide a 
final stage of matter-dominated-like oscillatory epoch, during which reheating can occur. A simple example for such 
a potential is 



V^Vo {bx^ + (x2)P) 



with 



& > 0, p> 1, 



(79) 



We have introduced this form for the potential because, for simple power-law monomials, i.e. V oc y^, with p > l, the 
second order action for the perturbations given in Eq. (1261) will vanish at y = since Q cx Vy = 0. This corresponds. 
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Figure 12. The evolution of 3MpiH^ /Vq for the potential V = Voie"^ — 1), and i/ = 1 (continuous black curve). This figure 
shows that after inflation ends (around A'^ ~ 72), the universe enters a matter dominated epoch, as the curve approaches a 
dashed line, which represents the line ln(3Mpi_H'^/Vb) — —3N + constant. This means that after inflation oc e"'^'^ oc a~^. 




Figure 13. Phase space plot for the variables dx/dt and x (left panel), and for dw/dt and w (right panel), during the oscillatory 
epoch. 

in general, to a strong coupling limit for the theory. One can avoid this situation by modifying the power-law potential 
as in Eq. ((79|) . There is no fixed-point at a; = for this form of the potential. Therefore, the field can oscillate around 
a; = to provide the mechanism to end the infiation without reaching Q = at y = 0. More in detail, according to 
the previous section, the points in region |x| > ^2/3 (unstable slow-roll fixed point of the dynamical equations of 
motion) will be forced to move to the region < •\/2/3, and the inflationary period will be long enough if the field 
X starts at = Xg > •\/2/3 with \w\ ~ 1, where w (x x' + "ix. The bottom line is that, in general, the 3-form field 
can drive long enough inflation without the ghosts or instabilities if its potential has local minimum at x = and has 
no local maximum between x — -iiXs- 

We also give another working example, the Gaussian potential, here defined as 

V = ^0(6"-^' - 1) , (80) 

which has similar properties to the power-law case discussed above. In fact, a long-enough slow-roll regime is followed 
by an oscillatory epoch where inflation ends. 
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Even if avoiding ghosts [Q > 0) and Laplacian instabilities (c^ < 0) are necessary conditions to be satisfied, they 
are not, however, sufficient, in general, to have a successful period of inflation. In other words, it is not assured that 
inflation ends for other classes of potentials which are, on the other end, free from instabilities. 

If the potential V{y) is such that for y > 0, it satisfies the conditions V > 0, Vy > 0, and Vyy > 0, then no 
instabilities arise, as already said. However, if we also impose that as y — >■ 0, we have V{y) — cy, where c is a positive 
constant, then for x ~ (and this point can be reached), Vcff oc x^, so that in general, an oscillatory epoch can take 
place, ending inflation. 

We have investigated the stability of the perturbations for a minimally coupled 3-form, whose action has a standard 
kinetic term and a generic potential function. We have found the conditions for which the inflationary dynamics 
can be stable, and gave some classes of potentials which can provide enough inflation without generating ghosts or 
Laplacian instabilities. We will leave the question to constrain the parameter space for this potentials by using the 
bounds on the spectral index and tensor-to-scalar ratio to a future research project. 
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Appendix A: The dual theory 



It is possible to define the 1-form or vector dual to the 3-form as 



A. 



(Al) 



where Eap-yS is the Levi-Civita antisymmetric tensor on curved backgrounds, which on Minkowski reduces to (.a/s-yS 
(with £0123 = 1 = -e^i^^). Then we also have E°'^^^ = e^'^^'^/v^. It is easy to show that \7 ^E^pyS = 0. In 
the following we will make use of the following relations e'^^'^^Cf^p^s = —6(5^, and e°'^'^^e^v^s = ^'2,{5'^5^ — ^"(5^). 
Therefore we obtain 



We also have 



48 



A^p^A^^-^ = e^p^se'-f-'^'B^B' = -dB^Bf^ . (A2) 
(ei23oVo-B° — eoi23V3i?'^ -I- e3oi2V2i3^ — e230iVii3^) 



X (ei230v"Bo - e°i23vi5j + £3012^2^^ _ ^2301^1^^) ^ ^ (V^B^)^ , 



so that the action is equivalent to the following one 



S 



(V^i?^)2 - V{Bl) 



(A3) 



(A4) 



which shows that the 3-form action is classically equivalent to a particular class of vector-tensor theories. Relation 
(|Aip can be inverted to give 



1 1 



07 ■ 



therefore, once the tensor A is known we can uniquely find B. At the level of the perturbations we find 



5B^' 



3! 7^ 



SA 



a/37 



(A5) 



(A6) 



which is valid on any background. 
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